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Friedel Oscillations as a Probe of Fermionic Quasiparticles 
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When immersed in a see of electrons, local impurities give rise to density modulations known as 
Friedel oscillations. In spite of the generality of this phenomenon, the exact shape of these modu¬ 
lations is usually computed only for non-interacting electrons with a quadratic dispersion relation. 

In actual materials, Friedel oscillations are a viable way to access the properties of electronic quasi¬ 
particles, including their dispersion relation, lifetime, and pairing. In this work we analyze the 
signatures of Friedel oscillations in STM and X-ray scattering experiments, focusing on the concrete 
example of cuprates superconductors. We identify signatures of Friedel oscillations seeded by impu¬ 
rities and vortexes, and explain experimental observations that have been previously attributed to 
a competing charge order. 
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I. INTRODUCTION 

In the study of strongly correlated materials, one 
of the common themes is the appearance of com¬ 
peting and/or intertwined orders (See for example 
Refs. [IHS] for an review). Recently this subject re¬ 
ceived considerable attention due to the experimen¬ 
tal observation of incommensurate charge modulations 
that coexist with high-temperature superconductivity in 
cuprateJ^i^. These modulations were originally ob¬ 
served on the surface of BSC CO via scanning tunnel¬ 
ing spectroscopy (STM]P^HIll, and recently found in the 


bulk of YBCO and other materials by X-ray scatter¬ 
ing experimentJ^HH, it is commonly believed that these 
modulations are adiabatically connected to the long- 
ranged charge order obse rved at large magnetic fields by 
quantum oscillationand nuclear magnetic resonance 
(NMR]PSI. In our earlier worlPH we proposed an alterna¬ 
tive interpretation: the modulations observed at small 
magnetic fields can be understood as Eriedel oscillations 
due to the scattering of quasiparticles with a short life¬ 
time, rather than as the evidence of a competing order. 

Eriedel oscillations are density modulations generated 
by local impurities acting on mobile charges, such as 
electrons in a metal. At the lowest order of perturba¬ 
tion theory, these modulations are proportional to the 
static density-density response function of the unper¬ 
turbed (homogeneous) system. Eor free electrons in three 
dimensions, this function can be analytically computed 
and its Eourier transform is peaked at twice the Eermi 
wavevector (see Ref. [22] for a review). As a consequence, 
Eriedel oscillations can b e exploited to directly measure 
the electron densit}®^^. In more complex materials the 
shape of Eriedel oscillations is determined by the band 
structure of fermionic quasiparticles, their lifetime, and 
the presence of a pairing gap. We suggest that the ob¬ 
servation of Eriedel oscillations is therefore a viable tool 
for studying the properties of quasiparticles in strongly- 
correlated materials. 

In this paper we theoretically analyze signatures of 
Eriedel oscillations in X-ray and STM experiments, focus¬ 
ing on the specific example of superconducting cuprates. 
The band structure of these materials has been ex¬ 
tensively studied by angle-resolved photoemission spec¬ 
troscopy (ARPES). In the normal phase, these materials 
have a single Eermi surface, whose phenomenological pa¬ 
rameters are well known. This information allows us to 
make quantitative predictions for the expected shape of 
the Ereidel oscillations. At low temperatures, the pres¬ 
ence of a superconducting gap challenges the main as¬ 
sumption of Eriedel oscillations, namely the presence of 
an underlying Eermi liquid. As we will see below, Eriedel 
oscillations can occur in a paired state as well, with some 
important differences. In particular, in a superconductor. 
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Friedel oscillations can be seeded by both local modula¬ 
tions of the chemical potential (as in normal metal) and 
local modulations of the pairing gap. 

In the following we fist present the theoretical frame¬ 
work used to analyze Friedel oscillations and then discuss 
its implications for recent experiments in cuprates. The 
appendices explain in detail several technical details of 
the calculations and in particular: how to p erfor m an 
ensemble average over static impurities (Sec. |A l| ); how 
to derive the formalism of r eson ant elastic scattering us¬ 
ing Green’s functions (Sec. A2) and h ow to relate it to 
the Lindhard susceptibility (Sec. |A3[ ); how to q uantify 
the effects of the spin-orbit co uplin g (Sec. A4) and of 
the light-held polarization (Sec. A5). 


II. THEORETICAL FRAMEWORK 


A. Lindhard susceptibility 


We open our theoretical description of Friedel oscil¬ 
lations by considering (hard) X-ray scattering experi¬ 
ments. This probe was used to detect Friedel oscilla¬ 
tions in vanadium-doped blue bronzd^, a charge density 
wave (CDW) material. For this material, accurate X-ray 
scattering experiments revealed two distinct incommen¬ 
surate diffraction peaks. These peaks where respectively 
identihed with the CDW wave-vector, and with Friedel 
oscillations at twice the Fermi wavevector. In a typi¬ 
cal X-ray experiment the intensity of the scattered light 
is proportional to the zero-frequency density-density re¬ 
sponse functiorP^ 

x(q)=^ dt J dx p{-x,t),p{-x,0) ), (1) 

where [•, •] is the commutation relation and p{x,t) = 
'0‘f‘(x, t)'0(x, t) is the charge density. For quasiparticles 
with a dispersion relation Sk and a hnite lifetime F, 
Eq. 0 becomes 




^k+q 


'^k '^k+q T 2ir 


( 2 ) 


where Uk = [l+exp((5/e—/i)/T)]“^ is the Fermi-Dirac dis¬ 
tribution function, and T the temperature. By neglecting 
interactions between quasiparticles, Eq. 0 disregards 
possible collective modes such as spin waves and para- 
magnons. In the case of free electrons (with 5k = /c^/2m 
and F ^ 0^ and at T = 0, Eq. Q can be evaluated 
analytically^. In two dimensions x(q) is momentum- 
independent for q < 2kF^ and decays algebraically for 
q > 2kF^ where kF is the Eermi momentunP^. In actual 
materials, the dispersion relation is more complex and an 
exact analytical evaluation of Eq. 0 is generically not 
possible. We therefore resort to a numerical evaluation 
of this expression. As we will see in Sec. Ill A[ this cal¬ 
culation leads to sharp peaks in x(g). When such peaks 
are observed in experiments, they may be interpreted as 
evidence of static charge density waves (CDW). 
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FIG. 1: Schematic diagram of a typical REXS experiment, in 
which X-ray photons scatter electrons from a core level to the 
conductions band (a), and viceversa (b). 


B. Local density of states 


In contrast to hard X-ray measurements, STM and 
REXS temporarily change the number of electrons in 
the conduction band and couple to the density of states, 
rather than to the density-density response function^. 
Specifically, at zero temperature the STM differential 
conductivity dI{T)/dV is proportional to the local den¬ 
sity of states ^(r, uj = V), given by the imaginary part of 
the retarded Green’s function: 


g{r,Lo) = Iin[G(r,a;)] 

■ pOO 

Jo 


( 3 ) 

( 4 ) 


Eor disordered materials, g{r,uj) varies in space and is in 
general unpredictable. It is therefore common to com¬ 
pute the two-dimensional Eourier transform of the signal 
at a fixed voltag^^ 

g{ci,L 0 ) = I d'^T g{r,L 0 ) . (5) 

As we will explain in detail below, the absolute value 
of ^(q, cj) depends on the types of scatterers present in 
the material, but not on their position (assuming that 
the sample is large enough to enable self-averaging of the 
scatterers’ position). 

Resonant elastic X-ray scattering (REXS) offers an 
alternative way to measure the local densiW of states, 
^(q, cj). As pointed out by Abbamonte et a/.“ STM and 
REXS describe analogous processes: in STM electrons 
tunnel to the sample’s conduction band from an atomic- 
size tip, while in REXS they are coherently pumped from 
a local core level (see Eig. [^. Based on this analogy, 
Abbamonted et al modeled the intensity of the REXS 
signal (at zero temperature) by 


^REXs(q, ^) 



duj'G^{uj — uj')g{q, uj') 


2 


(6) 


Here G^{uj) = [(cj + iFc)]”^ is the retarded Green’s 
function of the core level and Fc is its lifetime. In Ap¬ 
pendix A 2| we provide a derivation of Eq. 0 based on the 
Keldysh Green’s function formalism, which allows to ex¬ 
tend this expression to finite temperatures. Eurthermore, 
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in Appendix |A 3| we show that in the limit of non reso¬ 
nant scattering (Pc ^ oo) from a Fermi sea, Eq. § re¬ 
duces to the Lindhard susceptibility (©• Notably, Eq. (§ 
neglects the effects of the core-hole potential on the evo¬ 
lution of the conduction band. This effect is fundamen¬ 
tal to understand resonant inelastic scattering (RIXS)I^ 
processes, but can probably be neglected in the case of 
REXS. 

The prefactor A in Eq. ^ describes the transition am¬ 
plitude for the excitation of a single core hole. As shown 
in Appendix |A 4[ this quantity does not depend on the 
details of the core orbital and, in particular, is unaffected 
by the spin-orbit coupling. In the absence of magnetic 
impurities, the dipole approximation results into 


A oc {d\ {fji ■ r) (t)* ■ r) \d);, 


( 7 ) 


where \d) denotes the orbital wavefunction of the elec¬ 
trons forming the conduction band, r is displacement vec¬ 
tor in this state, and is the polarization of the incom - 


ing/outgoing photon. Eq. Q is used in Appendix A 5 


to 

compare the theoretical predictions of the present single¬ 
band model with the experimental results of Ref. [29] . 


C. Wannier functions and Bragg peaks 

We now consider the effects of non trivial Wan¬ 
nier functions on the Eourier-transformed local den¬ 
sity of states ^(q, cj). To achieve this goal, we first 
express Eq. in terms of the Eourier-transformed 
fermionic operators '0/c(t) = J d^r t) and 

their retarded Greens function G(k, k -|- q, t) = 
/y dt e-'‘^*([V>(k,t),?/)1'(k + q,t)]) as 

fl(q, w) = Im y] G'(k, k + q, w) , (8) 

k 

To derive Eq. ^ we assumed the system to be sym¬ 
metric under r — r. This symmetry allowed us to 
invert the order of the f d^k and Im operators (see SI-2 
of Ref. [30] ). This assumption is valid for example in the 
presence of a single scatterer at the origin of the axis. 
In the presence of several scatterers at random locations 
the present analysis applies to the abso lute value of the 
measured quantity (see Appendix |A 1| ). 

Eor a single-band model, the operator '0(r, t) is related 
to the annihilation of an electron (quasiparticle) on a 
single site, q, through the Wannier function IF(r — r^), 

=Y^W{r-ri)ci . (9) 

Combinin g Eq s. ^ and ^ we arrive to the 
expression®^ 

g{q, w) = Im y] WkG'iattice(k, k + q, w) Ilk+q , (10) 

k 

where W{]s) = f d‘^re^^'''W(r), G'iattice(k, k + q, i) = 

/y dt e“®“*([c(k,t),cl(k + qt)]), and Ck = 


giq.oi) 


1 
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FIG. 2: Schematic representation of the diagrams considered 
in the present analysis. Our approach is based on the first- 
order perturbation theory in the strength of the disorder and 
does not include the effects of interactions among quasipar¬ 
ticles. The function g{q,uj) is the Fourier-transformed local 
density of states. 


Note that by definition, Giattice(k, k-|-q, cu) is a peri¬ 
odic function of k and q with a period given by the 
Bravais lattice vectors, G. In the absence of impuri¬ 
ties Giattice(k,k + q,cj) = G(k,u;) X]Q(5(q - G), where 
G(k, cu) = G(k, k, cu). This expression gives rise to well- 
defined Bragg peaks in ^(q, cu), whose intensity is deter¬ 
mined by the width of the Wannier function. 

Cuprates posses a non-trivial Wannier function with d- 
wave symmetr 3 @. As explained in Ref. [35], this shape 
affects the phase and intensity of the Fourier-transformed 
STM signal at large wavevector (beyond the central Bril- 
louin zone), leading to subtle correlations that were in¬ 
terpreted as evidence of a competing order with d-form 
factoJ®^. In this paper we focus on the central Bril- 
louin zone, where the precise shape of the Wannier func¬ 
tion is not very important. For convenience, we then ap¬ 
proximate W (k) as a Gaussian wavefunction with width 
(7k = 1.8(27r/a), where a is the lattice constant. 


D. Scattering from local impurities 

In actual materials, as a consequence of disorder, the 
Fourier-transformed density of states ^(q, cu) is non zero 
even for wave-vectors that do not correspond to a lattice 
vector. Performing a first-order perturbation theory in 
the scattering potential (Born approximation) one findJ^ 

G(k,k + q,cj) =G(k,cj) ^5(q-G) 

G 

+ G(k, cj)T(k, q)G(k + q, cj) , (11) 

where G(k,co’) = G(k, k, cu), and T(k, q) describes the 
scattering of quasiparticles from momentum k to mo¬ 
mentum k + q. 

One of the main goals of this paper is to consider the 
effects of different types of impurities, defined through 
their scattering matrices T. We consider here only per¬ 
turbations that are static and quadratic in the quasi¬ 
particles’ creation and anihilation operators. Any such 
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perturbation can be described by the Hamiltonian 


Hpert — ^ ^ Q)^k^k+q (1^) 

k,q 


where T(k, q) = + 

giq-Xigik-(xj-Xi)^^ 

If the scatterer acts on a single site (or on the bonds 
linked to a single site), the scattering amplitude is given 
by the sum of two terms, which depend respectively on 
the momentum of the incoming and outgoing quasiparti¬ 
cles only: T(k, q) = Tk + ^k+q • Combining this expres¬ 
sion with Eq. we find 


5(q, a;) = ^ Im [WCG{k, a;) Wu+J ^ 5(q - G) 
k G 

+ ^ Im[Wk*G(k, uj){Tu + Tk+q)G(k + q, lo)W^+^ 


(13) 


Eq. (13) is at the basis of the present analysis: the first 
line corresponds to the density of states in an ideal lattice, 
while the second line describes the effects of the impu¬ 
rities (see also Eig. [^. In what follows we will mainly 
consider this latter contribution. 


E. Impurities in a paired state 


The above-mentioned formalism can be easily extended 
to include the effects of a spectral gap. Eor concreteness, 
we describe underdoped cuprates in terms of a single 
spectral gap, the paring gap A. Eollowing Ref. [30], we 
propose that the second energy scale observed in many 
experiments corresponds to the quasiparticle’s lifetime E, 
rather than to a distinct (competing) gap. This would ex¬ 
plain the uncertainty in determining the precise value of 
the gap in underdoped samples, varying roughly between 
A —T and A + T, depending on the type of experiment)^. 
The relatively-large value of I® in underdoped cuparates 
might be related to enhanced phase fluctuations, which 
lead to a loss of global phase coherence at the critical 
temperature Tc (see for example Refs. [39h47] b Notably, 
the density of states is a gauge-invariant object and, as 
such, depends only on the amplitude of A but not on 
its phase. Eor simplicity, we assume the pairing gap to 
have a pure d-wave form, Ak = Ao/2(cos/ca, — cos ky)^ 
although this assumption has little effect on the final re¬ 
sult. Conforming to the Born approximation, we assume 
Aq to be homogeneous over the sample and independent 
on the impurities. 

In the presence of a pairing gap, quasiparticles are con¬ 
veniently represented as 2 x 2 matrices in Nambu space 
(whose two entries are respectively particles and holes). 
In this notation the retarded Green’s function of a quasi¬ 
particle with momentum k and energy uj is given by 


G ^(k,u;) 


uj — 5k T /i T iT 

A_k 


Ak A 

U) -|- 5_k — /i + iT J 

(14) 


The dispersion relation 5k, the pairing gap Aq, and 
the quasiparticles lifetime T relevant to superconducting 
cuprates are provided in Sec. |m]and in Table [jl 


Static and quadratic perturbations can be divided 
into two main categories, CDW and PW, depending on 
whether they conserve the total number of quasiparticles 
(^ ctcj, see Eq. Jl2|)), or not (^ CiCj + cjcj), often re¬ 
ferred to respectively as diagonal and off-diagonal. In this 
paper we restrict our analysis to three specific types of 
impurities: two of them have a simple physical interpre¬ 
tation and correspond to local modulations of the chemi¬ 
cal potential (sCDW) and of the pairing gap (dPW). The 
third type (dCDW) corresponds to a local modulation of 
the intra-unit-cell nematic ordei® and has d-wave sym¬ 
metry. These three types of impurities correspond to the 
real-space Hamiltonians qCo,o, = 


(4,oCi,o+Co,oC-i.o)-(co,oCo,i+4.oCo-i)+H.c., = 


('^o.o'^i.o “I" ^o,o^-i,o) ('' 0 ,o'' 0 ,i “I" ^o,o^o,-i) “I" H.C.. The 

associated scattering matrices to be used in Eq. (13) are 


rpsCDW 

^k 



rpdCDW 

^k 


rpdPW 


dk 

0 

0 

dk 



(15) 


where dk = cos kx — cos ky . The Eourier transformed den¬ 
sity of states ^(q, cj) is obtained by numerically integrat¬ 
ing Eq. (13), with G(k, cj) and Tk respectively defined 
by Eqs. (14) and As we will see, a comparison be¬ 
tween the resulting plots and the experimental findings 
suggests a coexistence of sCDW and dPW, but rules out 
the presence of dCDW local modulations. 


III. EXPERIMENTS ON CUPRATES 

A. X-ray: effects of the band structure 


As mentioned in Sec. |H A , X-ray experiments cou¬ 
ple to the density-density response function and can be 
used to directly measure Eriedel oscillations. To em¬ 
ploy the Lindhard formula ([^ it is necessary to know 
the dispersion relation 5k, the chemical potential /i, and 
the quasiparticles’ lifetime T. In the case of super¬ 
conducting cuprates, these parameters can be directly 
read from accurate angle-resolved photoemission spec¬ 
troscopy (ARPES) experiments. Eollowing the common 
approach, we assume electrons to move within isolated 
CuO planes and map the conduction band in terms of 
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the two-dimensional dispersion relation 


Sk = ^ (cos kx + cos ky) cos kx cos ky 

+ ^ (cos kx + cos ky) (cos 2kx cos ky + cos kx cos 2ky) 

+ ^4 cos 2kx cos 2ky + ^ (cos 2kx cos kx + cos 2ky cos ky) . 

(16) 


In this work we specifically refer to three dis¬ 
tinct compounds: Bi2Sr2a:Laa^Cu06+^ (Bi2201), 

Bi2Sr2CaCu208+<5 (Bi2212), YBa2Cu307-x (Y123), 
whose band structure were experimentally deter¬ 
mined b y N orman et alP, Schabel et aim, Pasani & 
Atkinson^, King et The relevant parameters 

to — to are reproduced in Table |Tj Note that Y123 
material has inequivalent bonding (B) and antibonding 
(A) bands: Table refers only to the former one. 

The chemical potential /i is uniquely determined by the 
charge doping through the Luttinger count. Following 
the common convention, we denote by p the density of 
additional holes with respect to half filling: 


p = 2x — 1 , 


where 




(17) 


and k runs over the Brillouin zone. In the case of Y123, 
we identify the nominal doping with the algebraic average 
of the doping in the bonding and antiboding bands, p = 
{pa +Ps)/ 2. The resulting Fermi surfaces are plotted in 
Fig. ||;a). As shown in the inset, the Fermi surfaces are 
not circular, and display a significant amount of nesting 
at the ant inode 

Using the phenomenological parameters listed in Ta¬ 
ble [T| we can directly evaluate the Lindhard suscepti¬ 
bility (©■ Fig.lfb) presents x(q) along the direction 
(g, 0), for the three different materials. In all three cases, 
we observe a pronounced peak at a wave-vector rang¬ 
ing between 0.2 and 0.3. The exact position of the peak 
depends on the choice of the chemical potential, and is 
roughly equal to the distance between two adjacent antin¬ 
odes. The width of the peak is of order 0.03 — 0.1, lead¬ 
ing to a correlation length of about 10 — 30 unit cells, or 
40 — 120A. Its value is ma inly determined by the amount 
of nesting at the antinodes^ Among the three materials 
considered here, the sharpest peak is predicted in Y123, 
where the amount of nesting is maximal. In contrast, the 
Fermi surfaces of Bi2212 and Bi2201 involve a lower level 
of nesting, resulting in broader peaks. This could explain 
why, so far, (non-resonant) hard X-ray emeriments have 
revealed Friedel oscillations in Y123 only^. 

The specific choice of the band structure determines 
the details of the predicted signal. In the case of Y123, 
Fig. [^b) compares the signal resulting from the band 
structure of Pasani and AtkinsorP^ (continuous blue 
curve) and of Shabel et al^ (dashed blue curve). As 
shown in the inset of Fig. |^a), the latter band structure 
predicts a larger amount of nesting, in agreement with 


Band Structure 

Bi2212 

m 

Bi2201 

m 

Y123(B) 

m ISO] 

P 

0.0234 

- 0.148 

- 0.03 - 0.1256 

to 

- 0.5951 

- 0.5280 

- 0.42 - 1.1259 

ti 

0.1636 

0.2438 

0.1163 0.5540 

t2 

- 0.0519 

- 0.0429 

- 0.0983 - 0.1774 

ts 

- 0.1117 

- 0.0281 

- 0.353 - 0.0701 

U 

0.0510 

- 0.0140 

0 0.1286 

U 

0 

0 

o 

1 

o 

Doping, p 

0.04 

0.11 

Pb=-0.04 

( pA + pn )/ 2 = 0.12 

Lifetime, F 

0.004 

0.020 

0.002 0.001 

Gap, Ao 

0.040 

0.080 

0.030 0.030 


TABLE I: Phenomenological band structures used in this pa¬ 
per. The parameters to — ts correspond to hopping terms in 
a tight-binding model and are defined in Eq. |16| With re¬ 
spect to the originally published band structures, the chemi¬ 
cal potential has been shifted to achieve the required doping 
p (through the Luttinger count). Additionally, the parame¬ 
ter ts has been added to the band structure of Y123(B) in 
order to cure a spurious back-banding of the band structure 
(See SI-7 of Ref. [30]). F is the quasiparticle lifetime, used in 
Eqs. § and ( |14| ), and Ao the zero-temperature pairing gap 
used in ( [T4| . Their value can be read from the voltage de- 
pendence of the Eourier-tra nsform ed STM signaP^^^^^^®^^, 
or from ARPES experimentP^^. 


the experiment by Okawa et alW^. When used to predict 
the intensity of the REXS signal, the band structure by 
Ref. [51] predicts a peak at wavevector q = 0.28 with 
width dq ~ 0.05, while the band structure by Ref. [50] 
predicts a peak at q = 0.3 with dq ~ 0.03. For compari¬ 
son, the experiment of Chang et al¥^ shows a peak with 
maximal intensity at q = 0.31 and width 5q ^ 0.03, and 
is found to be in quantitative agreement with the present 
calculations. 

For completeness we mention that the peak observed 
in X-ray scattering experiments could be additionally en¬ 
hanced by effected that are not included in the present 
analysis. In particular, the Linhard formula <§ disre¬ 
gards the effects of the electron-phonon coupling. This 
coupling was instead found to be relatively strong in 
Y123 at this wavevector, leading to a significant phonon 
softenin^^. Electron-phonon coupling will generically 
lead to a sharpening of the X-ray response function, as 
well as to a renormalization of the position of maximal 
intensity. It seems plausible that the combination of the 
band structure of Ref. m with electron-phonon coupling 
could deliver a quantitative agreement with the experi¬ 
ments as good as the one obtained from the band struc¬ 
ture of Ref. [50] . As a side remark, we also note that the 
two predicted Lindhard response for Y123 differ by an 
overall multiplicative factor (see continuous and dashed 
blue curves in Fig. |^b)). This difference can be traced 
back to the different bandwidth predicted by the two 
models (^ 0.3eU in Ref. [51] and ~ leU in Ref. [50]) ]. 















6 



FIG. 3: (a) Fermi surfaces resulting from the band structures 
listed in Table(b) Lindhard response, Eq. along the line 
q = (27r/a) x (g, 0) for the same materials. 


Current experiments involve an unknown normalization 
factor and are therefore not sufficient to measure the ac¬ 
tual value of x{q) distinguish between these two sce¬ 
narios. Different experiments, and in particular resonant 
inelastic scattering (RIXS), might be able to fill in this 
information (see for example Ref. EZD- 

We now consider the full Lindhard susceptibility as a 
function of the two dimensional wavevector q = (27r/a) x 
{qx,qy)- Fig. El represents the results for Bi2212 (whose 
band structure is known to the highest degree of preci¬ 
sion) and displays three inequivalent local maxima. The 
global maximum occurs around the wavector = 

(27r/a) X (±0.5, ±0.5). This peak occurs in the other two 
materials as well (not shown) and its position is found to 
be independent on the doping level. Interestingly, 
corresponds to the wavevector of the anti-ferromagnetic 
order observed in the parent compound. Because the 





FIG. 4: Lindhard response, Eq.([^, as a function of q = 
(27r/a) X (qx^Qy) for Bi2212. The response of the other ma¬ 
terials is qualitatively similar (although the peak position is 
shifted away from q = 0.25). 


Lindhard formula describes both spin and charge sus¬ 
ceptibility, the predicted scattering enhancement around 
is a precursor of the long-ranged spin order achieved 
in the absence of doping (Mott insulator). 

A second broad peak appears at q = (27r/a) x 
(±0.25, ±0.25). The exact position of this peak is 
material-dependent and ranges between \qx\ = \qy \ = 0.2 
and \qx\ = \qy\ = 0.3 depending on the details of the 
band structure and the doping level, in analogy to the 
qy = 0 cut shown in Fig. I^b). Notably, this peak might 
easily escape experimentm probes: due to its broadness, 
it seems to merge with the stronger peak at q 7 r, 7 r, espe¬ 
cially if only the cut along the line q^ = qy is available. 
We will come back to this point in Sec. Finally, the 

third local maximum occurs at q = (±0.25, ±0.1): the 
wavevector q = (27r/a) x (±0.25,0) is predicted to be 
a saddle point sitting between these local maxima. We 
note that a similar behavior was observed in a recent ex¬ 
periments by Thampy et who found sharp peaks at 
q = (27r/a) x (0.25, ±0.015), separated by a saddle point 
at q= (27r/a) x (0.25, 0)P. 
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FIG. 5: STM spectra g(q, uj) along the line q = (^x, 0) x 27r for 
Bi2212 (see Table and T = j^dPW ^ ^ quasiparticle’s 
lifetime ((a) F = ImeV) leads to dispersive peaks, while a 
short lifetime ((b) F = 20meV) leads to non-dispersive peaks. 
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(a) SCDW (b) dCDW (c) dPW (d) sCDW+dPW (e) Experiment 
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FIG. 6: (a-d) Numerical evaluation of Fourier-transformed STM spectra, Eq. (13), for different types of impurities: local 
modulation of the chemical potential (sCDW), local modulations of the intra-unit-cell nematic order (dCDW), local modulations 
of the pairing gap (dPW). Model details: see column Bi2201 of Table [T| Th e black circles denote the Bragg peaks at q = (±27r, 0) 
and (0,±27r) . (e) Experimental measurement reproduced from Ref. [29]. 


(a) sCDW (b) dCDW 


(c) dPW 


(d) sCDW+dPW 


(e) Experiment 




H (Inla) 


EIG. 7: (a-d) Numerical evaluation of the REXS signal, Eq. for different types of impurities (see caption of Eig. [^. In 
order to easy the comparison with the experimental plot, the theoretical predictions in the (g, g) direction (orange diamonds) 
have been reduced by a factor of 0.25 with respect to the (^,0) direction (green squares). Model details: see column Bi2201 
of Table and Fc = 300meV. (e) Experimental measurement reproduced from Ref. m- The experiment ally-observed peak at 
q ^ 0.25 in the (^,0) direction (and its absence in the {q,q) direction) is correctly reproduced by our simple model of Eriedel 
oscillations. 


B. STM: dispersive vs non-dispersive peaks 


We now proceed to discuss STM experiments, by first 
offering a brief summary of the main results of Ref. m- 
Specifically, in that paper we related the emergence of 
non-dispersive peaks in underdoped cuprates to their rel¬ 
atively large inverse quasiparticle lifetime F. In materi¬ 
als where F is small (such as overdoped cuprates), the 
STM probe excites quasiparticles with an energy that 
precisely corresponds to the tip-sample voltage. In this 
case, energy and momentum conservation leads to the 
well-known “octet model”!^. This model predicts the 
emergence of seven inequivalent dispersive peaks, which 
can be found by connecting points on the Fermi surface 
where the pairing gap is equivalent to the tip-sample volt¬ 
age. As shown for example by Nowadnick et these 
peaks are indeed reproduced by Eq. (13) in the limit of 
F 0. In contrast, for a finite F, the argument leading 
to the octet model does not apply because the quasipar¬ 


ticles’ energy is not conserved. In this case, a numerical 
evaluation of Eq. (13) is necessary. As shown in Ref. m 
these calculations lead to non-dispersive peaks around 
the wavevectors connecting the antinodes. These scat¬ 
tering wavevectors are enhanced at all voltages for two 
reasons: (i) Any scattering is enhanced at the antinodes 
due to the Ferm i surface nesting (in analogy to the analy¬ 
sis of Sec. Ill A[); (ii) The modulations of the pairing gap, 
j^dPW {Ib)^ are proportional to the pairing gap Ak 

and are therefore enhanced at the antinodes, where the 
latter is maximal. 


The effect of F on the calculated STM maps is high¬ 
lighted in Fig. where F varies from ImeV (subplot 
(a)) to 20meV (subplot (b)), while all other parameters 
are kept fixed. The former plot displays dispersive peaks, 
while the latter mainly non-dispersive ones. Importantly, 
the temperature dependence of F can explain the transi¬ 
tion between dispersive peaks (at low temperatures) and 
non-dispersive peaks (at higher temperatures) reported 
in Ref. [60]. 














































C. STM and X-ray: identifying the impurities 

To further clarify the nature of the main source of 
disorder (sCDW,dCDW, or dPW), we now consider the 
STM and REXS experiments of Comin et a/P. Their 
two-dimensional Fourier-transformed STM signal is re¬ 
produced in Fig.j^e). The intensity of the signal is max¬ 
imal in a cross-shaped region, oriented in the (±g, 0) and 
(0, :hq) directions. Figs, [^a-c) represent our theoretical 
calculations for the three types of impurities defined in 
Eq. ( p!^ . The correct shape of the signal is reproduced 
only by local modulations of the pairing gap (dPW - sub¬ 
plot (b)), suggesting that this is the dominant sources of 
disorder. A similar conclusion was reached by the inde¬ 
pendent analysis of Nunner et al^. The experimental 
REXS measurement of the same material is reproduced 
in Fig. [^e). It shows a pronounced peak in the (^,0) 
direction, and a monotonous behavior in the {q, q) direc¬ 
tion. A comparison with the theoretical curves. Fig. [^a- 
c), reveals that this effect is reproduced only by local 
modulations of the chemical potential (sCDW - subplot 
(a)). 

This analysis leads to an apparent inconsistency: STM 
reveals local modulations of the d-wave pairing gap 
(dPW), while REXS reveals local modulations of the 
chemical potential (sCDW). The solution of this appar¬ 
ent paradox is hidden in the intrinsic properties of the 
two probes: STM measurements refer to low voltages and 
probe the scattering of quasiparticles with small energy 
E ^ Ao ~ 20meV. In contrast, REXS probes the scat¬ 
tering of quasiparticles with energy F < Pc ~ 300meV. 


Due to the coherence factors appearing in Eq. (13), 
quasiparticles at different energies are mainly affected 
by different sources of disorder: low-energy quasiparti¬ 
cles are mainly affected by modulations of the pairing 
gap, while high-energy quasiparticles are mainly affected 
by modulations of the chemical potential (see also SI-3 
of Ref. [30] ). This effect becomes evident in the present 
calculation: the intensity of the STM signal is signif¬ 
icantly stronger for dPW (Fig. [^c)) than for sCDW 
(Fig.m: a)), while the intensity of REXS is stronger for 
sCDW (Fig. [^a)) than for dPW (Fig. [^c)). The ex¬ 
perimental results are then best reproduced by a super- 

g )sition of both types of modulations (Figs. I^d) and 
d)). This result is in line with Jeljkovic et alP^, who 
found a strong correlation between local perturbations 
of the pairing gap and of the chemical potential (identi¬ 
fied there as atypical oxygen vacancies). Notably, local 
modulations of the intra unit-cell nematic order (dCDW) 
are inconsistent with the q dependence of both STM and 
REXS signals. 

Let us now discuss a theoretical prediction made in 
Ref. [30], which appears to be in contradiction by the 
experiment of Comin al^. Specifically, Ref. ISO] pre¬ 
dicted the existence of a peak in the REXS signal at 
wavevector (0.25,0.25). In contrast, the experimental 
measurements of Ref. [29] (orange curve in Fig.j^e)) does 
not show any significant peak at (0.25,0.25). We believe 





FIG. 8: Two dimensional plot of the predicted REXS sig¬ 
nal. Pronounced peaks are observed at q = (27r/a) x (0.25, 0) 
and q = (27r/a) x (0.5,0.5). Model details: Bi2201 with 
sCDW+dPW impurities. 


that the absence of the peak at (0.25,0.25) is due to its 
blending with the larger and broader peak at (0.5, 0.5) 
along the same direction.!^ This phenomenon is clearly 
demonstrated in Fig.[^ showing the predicted REXS in¬ 
tensity as a function of the two-dimensional wavevector 
q. Note the close analogy with the results of the Lind- 
hard susceptibility shown in Fig.|^ We hope that future 
experiments will be able to confirm the hereby prediction 
of an increased scattering in the (q,q) direction. 


D. X-ray: c-axis correlations 

Until this point we considered two-dimensional models 
only and analyzed correlations along the a and b principal 
directions only. Recently, Gerber et found that the 
c-axis correlations provide a clear distinction between the 
short-ranged modulations observed at low magnetic fields 
and the long-range modulations found at large magnetic 
fields. Specifically, while the former is peaked around 
kz = 0.5 X (211 jc) (or equivalently has a period of two 
unit cells), the latter is peaked at integer wave-vectors. 
As we will now explain, this observation is consistent 
with Friedel oscillations seeded by an impurity sitting 
at the interface between two unit cells. This situation is 
naturally realized in Y123, where the CuO chains are the 
main source of disorder and are equally spaces from the 
two neighboring CuO planes (see Ref. [65] for a review). 

If we neglect the tunnelling of electrons in the c direc¬ 
tion (i.e. among planes belonging to distinct unit cells), 
we obtain electronic bands that do not disperse in this 
direction. If we additionally assume that the scatter¬ 
ing matrix T and the Wannier function W are separable 
functions of the spatial coordinates, the qz dependence 
can be factored out from Eqs. @ , ( [T^ and (§, leading 
to /REXs(q) = where 

9z{qz) = Y.WlT{K,q,)Wk^+,^ (18) 

k. 
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(a) Symmetric defect 


Anti-symmetric defect 



CuO chain 

CUO 2 plane 
CUO 2 plane 

CuO chain 
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FIG. 9: Crystal structure of two neighboring unit cells of 
YBCO. A symmetric (anti-symmetric) defect located at a 
CuO chain induces the same (opposite) charge displacement 
on the CuO 2 planes of the two unit cells. Upward (down¬ 
ward) arrows indicate increased (decreased) charge density, 
graphical representation from Ref. |63l. 


(locally) on two neighboring CuO planes only: 


5i^o ± 6i^i 


(19) 


Here + (—) refers to a symmetric (antisymmetric) im¬ 
purity, the index i runs over the CuO layers, and c is 
the lattice vector in the z direction. A more accurate 
description should take into account the distinction be¬ 
tween bonding and anti-bonding bands, but we defer this 
point to a future study. The effects of symmetric and 
anti-symmetric impurities are schematically plotted in 
Fig.U The Fourier-transformed scattering amplitudes 
of symmetric and anti-symmetric impurities are then re¬ 
spectively given by 


T^{kz,qz) = 2cos 1 

(?)• 

(20) 

T~{k^,q^) = 2sin ( 

:f)- 

(21) 


Note that for symmetric impurities the intensity of T 
is peaked at integer kz{c/27T) = 0,±1,±2,.., while for 
antisymmetric impurities it is peaked at half integer 
kz{c/27r) = ±0.5, ±1.5,.... The intensities of the X-ray 
signals at low and high magnetic fields are therefore re¬ 
spectively consistent with anti-symmetric and symmet¬ 
ric impurities. To reproduce the experimental observa¬ 
tions, we need to introduce trial Wannier functions. For 
simplicity we again refer to Gaussian function W{z) = 
e~^ /2c which allow an analytic evaluations of Eq. 
leading to 
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q, (c/27r) 



q, (c/27r) 


and 


9z{qz) 

= (f) 

(22) 

9z{qz) 


(23) 


respectively for symmetric and antisymmetric impurities. 
These curves are plotted in Fig. 10 subplot (b) repro¬ 
duces the position and width of the experimentally ob¬ 
served signaP^. 


FIG. 10: Intensity of X-ray scattering as a function of the 
c-axis component of the wavevector Qz for (a) symmetric and 
(b) anti-symmetric impurities, Eqs. ( |23| ) and (23). The latter 
curve is similar to the signal observed by Gerber et a/.® at 
small magnetic fields and peaked around Qz = ±1/2. 


and iREXsjkx, ky^uj) has been computed in the previous 
sections. Eq. (1^ demonstrates that a non-trivial qz de¬ 
pendence can be obtained due to the shapes of the impu¬ 
rity and of the Wannier functions of the conduction-band 
electrons. The same approach goes through for hard X- 
ray diffraction, where one finds x(q) = g{qz)x(kx^ ky^uj)^ 

• i 


with x{kx,ky,uj) given by Eq 
pendix A3). 


(see Eq. A13 in Ap- 


Eor simplicity we now focus on an impurities that act 


IV. SUMMARY AND OUTLOOK 

In this paper we presented a theoretical modeling of 
recent X-ray, REXS, and STM measurements of under¬ 
doped cuprates, with specific attention to Ghiringhelli et 
a/.® Chang et alW, Comin et al^. To interpret their ex¬ 
perimental findings, these authors assumed the existence 
of a competing order, distinct from superconductivity, 
and associated with the spontaneous breaking of trans¬ 
lational symmetry. The pseudogap energy scale could 
then correspond to the excitation gap required to re¬ 
store the translational invariance. The association be¬ 
tween the charge ordering and the pseudogap phase is 
however undermined by recent X-ray experiments reveal¬ 
ing the same type of charge ordering in electron-doped 
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cuprate^, where a pseudogap phase is not expected to 
subsist. In addition, the experiment by Gerber et 
showed that the oscillations observed at small magnetic 
fields have different c axis correlations than those ob¬ 
served at large magnetic field^, indicating that these are 
two distinct effects. Similarly, recent measurements of 
the Hall conduct ivit}EI showing that the long-range or¬ 
dered modulations appear only for magnetic fields that 
are larger than a critical value He ~ 20T. 

Following Ref. [30] we propose here that the modula¬ 
tions observed at small magnetic fields are simply due 
to Friedel oscillations around local sources of disorder. 
Because our interpretation is based on the Born approx¬ 
imation (first-order perturbation theory in the impurity 
strength) and we consider each scatterer independently, 
we expect the correlation length of the modulations to 
be independent on the concentration of impurities. This 
prediction has been now confirmed by two experimen¬ 
tal observations: (i) Achkar et modified the amount 
of disorder in Y123 through a thermal quench and ob¬ 
served that the correlation length of the observed oscil¬ 
lations was unchanged, (ii) The analysis of materials 
with similar band structure and different amount of in¬ 
trinsic disorder (such as Bi220J^, Bi220J^, Bi221^, and 
Hgl20 revealed an approximately constant correlation 
length. These findings are not consistent with theories 
of competing orders, in which the predicted correlation 
length should be directly related to the amount of exter¬ 
nal disordeJ^. 

By considering the scattering of short-lived quasipar- 
tice from local impurities, we can quantitatively repro¬ 
duce all the experimental findings: Our model correctly 
predicts the wavevector and correlation length of the spa¬ 
tial modulations that were observed in X-ray (Figs, [^b) 
and|^, STM (Fig. [^, and REXS (Fig. experiments. 
The wavevector is similar (but not identical) to the dis¬ 
tance between adjacent antinodes, where the Fermi sur¬ 
face is often quite nested. Our approach reproduces ex¬ 
perimental observations that were interpreted as evidence 
for the d-wave symmetry of the oscillations (Figs.[^- see 
also Ref. [35] for an in-depth analysis of the phase corre¬ 
lations observed by Fujita et al)^. Finally, it naturally 
accounts for the non-trivial c-axis correlations observed 
in X-ray experiments (Fig. [Tq|). 

To reproduce the experimental results, we introduced 
different models of local impurities and found that the 
most dominant type corresponds to local modulations of 
the chemical potential and of the pairing gap. In STM 
maps, the former contribution is generically dominant 
along the (q,q) direction, while the latter is dominant 
along the (0,q) direction (see FigJ^. The interplay be¬ 
tween these two sources of disorder connects to the earlier 
analysis of STM data in the presence of a magnetic field 
performed by Hanaguri et al^ and He et alP^. These 
authors found that the ratio between the (q,q) and ( Q,q ) 
components generically increases with magnetic fielcP^. 

This effect can be understood by noting that in type- 
H superconductors external magnetic fields generate iso- 


T 



FIG. 11: Proposed phase diagram of BSCCO and YBCO 
compounds, following Ref. m-- the long-range-ordered CDW 
phase observed at high magnetic field is distinct from the 
pseudogap (PG) phase, characterized by an incoherent pair¬ 
ing gap of preformed pairs. The inset refers to doping levels 
p 0.1, where a direct transition between the superconduct¬ 
ing (SC) and CDW phases is observeeP^. 


lated vort exes, in whose core the pairing gap is locally 
suppresseeP^^^®^. Vortexes are then similar to other 
types of local impurities, and generate Friedel oscillations 
around their center. This effect was for example directly 
observed by Simonucci et who found Friedel os¬ 

cillations around magnetic vortexes in the self-consistent 
solution of the BCS equations in the presence of a vortex. 

When the density of vortexes reaches a critical value, 
they can depin from local defects and give rise to a long- 
ra nge-o rdered phase. Following the proposal of Wu et 
a/j 2Q | 75 ( we believe this effect to be responsible for the 
formation of a long-range-ordered phase at large mag¬ 
netic fields that was observed by quantum oscillations^, 
NmfP^RH^ and sound velocit}C3 experiments. Indeed, 
the measured critical field ^ 20T corresp onds t o an av¬ 
erage distance between vortices of d = ^ lOOA, 

which is comparable with the correlation length of Friedel 
oscillations. The CDW phase observed in cuprates would 
then be analogous to the field-induced spin density waves 
(FISDW) observed for example in Bechgaard salts (see 
Ref. [78] for a review). Because the correlation length 
of Friedel oscillations depends on the amount of nesting 
at the antinodes, it is natural to expect the magnetic 
phase to be enhanced around p = 0.1, where the antin¬ 
odes are maximally nested. The resulting phase diagram 
is plotted in Fig. IE and highlights our claim that the 
long-range ordered CDW phase is distinct from the pseu¬ 
dogap (PG) phase observed at zero magnetic field. 

We now discuss how to utilize STM maps to further 
compare the effects of magnetic fields, temperature, and 
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(a) small H 


(c) small T 


(e) large p 
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g(q, 3mV) 


V V = 16meV T = 4.5K -f 


(b) large H 


(d) large T 


(f) small p 



(g) dominant direction 



( 1 . 1 ) 


(- 1 . 1 ) 
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FIG. 12: Fourier-transformed STM measurements of different materials: (a-b) Ca2xNaccCu02Cl2 (Tc=28K) at low and high 
magnetic field. Reproduced from Ref.^ (c-d) Bi2212 (Tc = 37K) at low and high temperature. Reproduced from Ref. [73] : 
(e-f) Pb-Bi2201 at large and small hole doping. Reproduced from Ref. [711 (d) Deep in the superconducting phase the main 
source of scattering is along the (1,1) and (1,-1) directions, (h) When approaching the pseudogap phase the scattering is mainly 
along the (0,1) and (1,0) directions, signaling the presence of phase inhomogeneities. 


doping. Fig. shows that approaching the pseduo- 
gap phase (by increasing the temperature, or decreasing 
the doping) generically leads to an increase of scatter¬ 
ing in the (0, q) direction, which is associated with local 
modulations of the pairing gap. This observation is in 
agreement with recent muon spin rotatiorP^ (/iSR) and 
NMRP^ experiments, which detected enhanced static in¬ 
homogeneities in the pseudogap phase. A similar conclu¬ 
sion was reached in Ref. [80], where the effects of disor¬ 
der were found to be similar to the effects of temperature 
and magnetic fields (see also Ref. [81], where a pseudo¬ 
gap phase was found in disordered thin films). Inhomo¬ 
geneities of the pairing gap are naturally accompanied by 
a reduction of the long-range coherence: the transition 
to the pseudogap phase may be due to a loss of coherence 
of the pairing gapF^, rather than to its disappearance. 

From the prospective of fermionic quasiparticles, the 
transition to the pseudogap phase is generically asso¬ 
ciated with an increase of the inverse lifetime F. The 
role of this quantity on ARPES measurement is well 
known and offers a simple explanati on for t he “ Fermi 
arcs” observed in underdoped cuprate^^^^^^^^. STIVP^^ 
and transporfpS measurements show that the inverse 
quasiparticle lifetime F is strongly enhanced in under¬ 
doped cuprates and probably diverges at the transition 
to the Mott insulator. This observation suggests a pos¬ 


sible relation between F and the critical temperature of 
cuprateP^. Measuring the temperature dependence of F 
would allow to distinguish between the effect of disorder 
(elastic scattering, which exists down to zero tempera¬ 
ture) from the effects of interactions (inelastic scatter¬ 
ing, which is supposed to diasppear at zero temperature). 
Surprisingly, although the inverse quasiparticle lifetime 
F is commonly used in fitting virtually any experimental 
spectroscopic data, a systematic study of this quantity 
as a function of doping, temperature, and magnetic field 
has not been performed yet. We hope that the present 
work will motivate a new analysis of existing data along 
these lines. 
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Appendix A: Technical details 
1. Average over impurities 


In this section we consider the effects of several identi¬ 
cal impurities, located at random positions. We find that 
the absolute value of ^(q, cj) is independent on their posi¬ 
tion and therefore is an intrinsic property of the system. 
In the text, Eq. we assumed g(r,cu) to be symmet¬ 
ric under r ^ — r. This approximation is valid in the 
presence of a single impurity located at the origin of the 
axi^SH Xo extend this treatment to systems with sev¬ 
eral impurities, we first notice that within the present 
framework (first order perturbation theory in the impu¬ 
rity strength), g(q ^ G,Cc;) is given by a sum of terms, 
each referring to the scattering of quasiparticles from a 
single impurity: ^(q, cj) = X^^^i(q, cj), where % runs over 
all the impurities. To compute ^^(q, cj) we first consider 
a coordinate system whose origin is located at the center 
of the ith impurity, where Eq. applies. We then shift 
back ^^(q, cj) to the common lab frame by the multipli¬ 
cation with and sum all the terms. In the case of 

N identical impurities we obtain 

N 

^(q, w) = E ■ (Ai) 

i=l 


Here go{q^uj) is the scattering amplitude from an iso¬ 
lated impurity locat ed a t the axis origin, computed from 
Eq. ( [l^ . In Eq. (AI) the phase of ^(q,cj) is deter¬ 
mined by the random positions and is therefore not 
predictable. In contrast, its absolute value averages to 


(l5(q,w)|) = Ngo{ci,Lo) . (A2) 

Here we used the observation that by definition, go{q,uj) 
is a real function. 


2. Green’s function approach to REXS 

In this section we derive Eq. (§ within the Keldysh 
path-integral formalism. This approach allows us to ex¬ 
tend the results of Abbamonte et alWi to finite temper¬ 
atures. In REXS experiments X-rays are scattered upon 
the material to be examined at a frequency that allows 
the creation of a core hole, i.e. the excitation of an inner 
orbital of the atom to the conduction band (see Fig. [^. 
The action of the incoming field can be described as 
Vin = Eine^^^6{t)d^H.C., where Ei^ > 0 describes the 
amplitude of the incoming x-rays, cj its frequency, d and 


c are fermionic operators describing electrons (quasipar¬ 
ticles) respectively in the conduction band (in cuprates 
formed by d orbitals) and in the core level. Shortly after, 
an electron from the conduction band fills in the core level 
and emits an X-ray photon, which is observed by the ex¬ 
perimental setup. This decay process can be described by 
the operator Vout{t) = d^ where creates 

an outgoing photon and A is the light-matter coupling. 
In perturbation theory, the outgoing field is given by 

Eont{t) = (aout) - (Jde-"'"* + H.c.) (A 3 ) 


where we assumed the initial state to be empty of outgo¬ 
ing photons. Applying perturbation theory (and neglect¬ 
ing oscillating terms), we obtain the Kubo formula 

Eout{t) = iEinQ{t){[d\o)c{0),c'^{t)d{t)]) (A4) 


where 0(t) is the Heaviside theta function and 
i s th e commutation relation. In Keldysh notation, Eq. 
(A4) becomes the sum of eight terms with an odd number 
of “classical” fields c, d and of “quantum” fields c, d. Four 
of these terms contain three quantum fields and their 
expectation values are identically equal to zero. We are 
then left only with terms containing three classical fields 
and one quantum field: 


Eout (t) =^Eine{t) (f ( 0 )c( 0 )c* {t)d{t) 

+ d*{Q)c{Q)c*{t)d{t) + (i*( 0 )c( 0 )c*(i)d(i) 

+ d*( 0 )c( 0 )c*(t)(l(t)) (A 5 ) 


For t > 0 only the first two terms are non zero (and for 
t < 0 the last two are non zero). Eq. ( |A5[ ) can be fur¬ 
ther simplified by introducing the retarded and Keldysh 
Green’s functions = (d*d) and = (d*d): 


= iEi„ [G^{t)G^it) + Gf (i)Gf (t)] (A 6 ) 


Each of the two terms of Eq. (A6) corresponds to the 


product of two Greens functions, evaluated at the same 
time, or equivalently their convolution in the frequency 
domain: 


• roc 

Eoutiuj) = ^Ein duj'[G^{uj-uj')G^{uj') 

(cj')Gf (cj - cj')] . (A7) 

At thermal equilibrium the Keldysh components sat¬ 
isfy the fluctuation-dissipation theorem G^(cj) = 
2Im[G^] tanh(cj/2T) « 2Im[G^]sign(cj) and G^ = 
2Im[Gf] tanh((cj —E/,,)/2T) —2Im[Gf]. In the limit of 

T ^ 0 we find 

/ OO 

dio' [lm[G^]{Lo-Lo')G^{Lo') 

-OO 

+ sign(:^')Im[Gf](c^')Gf - oj')] . (A8) 
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The real component of Eout (the component that is in 
phase with Ein) has a particularly simple form 

poo 

Re[K«t](w) = 2Ein / dLo'lm[G^]{LO - Lo')lm[G^{Lo')] 
Jo 

(A9) 

Applying Karmers-Kronig relation we then obtain 

poo 

Eoutito) = 2Ein / dco'G^ico - Lo')Im[G^{Lo')] (AlO) 
Jo 

For a featureless core level with response function 
G^{q,uj) = [(cj + iFc)]”^, we recover exactly the same 
expression as in Ref. [27] and Eq. with A = 2£’in. 


3. Prom REXS to Lindhard 


In the main text we provided an expression for the in¬ 
tensity of the REXS signal at zero temperature, Eq. Ill- 
Here we show that, in the case of non-resonant scat¬ 
tering from a Eermi gas this expression simply reduces 
to the Lindhard susceptibility The present deriva¬ 
tion is a corollary of a more generic relation between the 
non-resonant limit of resonant inelastic scattering (RIXS) 
and density-density response functions (see Ref. [85] for 
a review), and is brought here for completeness. Non¬ 
resonant scattering can be described as a REXS process 
in the limit of Tc ^ oc. In a Eermi gas with a local 
impurity, Go{k,uj) = l/{uj — Sk E '^^0+), W{k) = 1, and 
T{k^k E q) = 1. Under these conditions Eq. in and 
give IreXS -fxray = \ {AC {(\) /Tc\^ With 


pOO 

C(q) = / dhj' y' Im 

Jo ,, 


T — 'S^k+q T 

(5(w' - ek+q) , d{uj' - ek) 


OJ' — £^k+q 


Jo E 

-E 


k 

77/k 


'^k+q '^k 

^k+q 


'^k '^k+q 


^k+q ^k ^k ^k+q 


(All) 


4. Spin-orbit effects in REXS 

In this section we study the dependence of REXS scat¬ 
tering on the polarization of the incoming (i) and outgo¬ 
ing (o) photons. As an important result, we will show 
that in the absence of magnetic impurities, the intensity 
of the REXS signal is not affected by spin-orbit effects. 
Eor an isolated atom, the REXS intensity / is given by 
the product of dipole matrix elements for the absorption 
and the emission: 

lifiufio) oc {fil ■ {xl)i\v\^n)) iVi ■ (i/>„|r|i/ii)), (A14) 

where Ei is the initial (and final) core electron state and 
jpn is a valence Mx 2 _y 2 orbital with spin a at the same 
site. In typical experiments one selects a resonance so 
that only 2p core levels with total angular momentum 
j = 3/2 are excited to the valence band. Thus the to¬ 
tal polarization-dependent intensity is the sum over spin- 
orbit eigenstates rrij = —3/2, —1/2,1/2, 3/2: 

oc 5](r):.(2p^/2|r|34^_,.,<T)) 

ruj 

X {fji ■ {3d^2_y2,a\r\2pf^E > (^15) 

To compute the dipole matrix elements, we introduce a 
unit operator in the basis of separate spin and orbital 
angular-momentum eigenstates 

/ oc E ( rj* ■ {2p^f\mi,ms){rni,ms\r\3dx2-y2,a 

rrij , 1711 ,iris ^'m'g,im'g 

X [rii ■ (342_y2,cr|r|m(,,m0(m(,,m;|2py2^) 

(A16) 

= E {nl-{2‘PE\rni,(T){rni\v\3d^2_y2)^ 

rrij 

X (rji ■ (3d^2_y2\r\m'^}(m'^,a\2pKE (^17) 

The Clebsch-Gordan matrix elements vanish unless rrij = 
mp E cr = m'p E cr^ and hence we require rrip = mp. We 
then have the further simplification: 


Here in the transition from the first to the second line we 
used l/{x E i0+) = 1/x — i7rS(x), and in the transition 
from the third to forth S(u; — Sk) = rik, where rik is 
the Eermi-Dirac distribution at T = 0. We obtain 


I oc ^{2pflf\me,a) 

rrij ,m£ 


( 77 * • {me\r\3d^2_y2)) 


X {fji ■ {3dx2_y2\r\me)) 


(A18) 


-^Xray 



^k ^k+q 
^k ^k+q 



(A12) 


The present derivation can be extended to the case of 
non-trivial Wannier functions, and scattering amplitudes 
of the form T(k, k + q) = Tq leading to 


-Ucray 


A 

f; 


^tRkW+q 


^k ^k+q 

'^k '^k+q 


2 


(A13) 


If the Hamiltonian is spin-independent, i.e. if there is no 
spin-density wave, the amplitude is independent of the 
spin a of the photoelectron in the intermediate state and 
thus the two spins contribute equally to the coherent sum 
over histories, and we have 


rrii rrij,a 

X {vi ■ {3d^2_y2\r\me)) 


{Vo ■ {me\r\3d^2_y2)) 


(A19) 
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{2VrL]\mt,a) 


2 

is the probability that a core 


electron with orbital angular momentum mg and un¬ 
known spin is in a total spin-j = 3/2 state. By spherical 
symmetry this is obviously independent of since mg 
is coordinate-dependent but j is not. Since this is an 
m^-independent quantity, we obtain 


/ oc E (C ■ W|r|342_y2)) {f]i ■ {3d^2_y2\r\me)) 

mi 

(A20) 

=( 342 _y 2 | {fji • r) {fjl • r) | 342 _j, 2 ) (A21) 

We now have a tensorial matrix element that is not mod¬ 
ulated by the spin-orbit effect except for the aforemen¬ 
tioned constant prefactor that represents the contribu¬ 
tion to resonant scattering only from j = 3/2 core states. 


5. Polarization dependence of REXS 




In this section we study the dependence of the REXS 
signal on the wavevector of the incoming photon k. This 
dependence was experimentally measured by Comin et 
a/.“ and used to identify the dominant type of charge 
modulations. According to the present single-band ap¬ 
proach, the predicted k dependence is instead identical 
for all types of modulations: Our model corresponds to 
the 5-wave case considered by Comin et al. . This model 
is found to be in good agreement with the experimental 
measurements (see Fig. [T^. 

Our starting point is Eq. 0 . Because the outgoing 
beam is not filtered according to its polarized, the mea¬ 
sured signal is proportional to the sum of the intensities 
of the two outgoing polarizations: 


^REXS 

o=a' ,7 t' 


(A22) 


where the tensor M is defined by = {d\ro^r^\d)^ 

i = cr, TT is the incoming polarization, and o = cr',7r' is 
the outgoing polarization. We denote by F the diago¬ 
nal matrix corresponding to M in the principal axis of 



FIG. 13: Experimental setup of Ref. [29]. 


FIG. 14: Polarization dependence - comparison between the¬ 
ory and experiment. The continuous curves correspond to 
Eq. (A23) and Eq. (A24), obtained from Eq. ( A22| in the 
case oi Fz =0. The dotted curves correspond to Eq. ( |A22 ) 
in the case of Fz = 0.15. The dashed curves are reproduced 
from Fig. S5 of Ref. [29] and include the corrections due to 
self-absorption. The errorbars indicate the experimental mea¬ 
surement of the Y675 sample reported in Ref. [29]. 


the lattice. Its three non-zero entries are Fx = (d|x^|d), 
Fy = {d\y‘^\d), and F^ = {d\z‘^\d). The ratio between 
these quantities was measured in Ref. m and found 
to be FzjFx ~ FzjFy ^ 0.15. The smallness of F^ 
indicates that the conduction band has a small exten¬ 
sion in the 2 : direction, in agreement with the theoretical 
calculation^, which predict a dominant charac¬ 

ter. In the experiment, the direction and the polariza¬ 
tion of the incoming photons are kept fixed, while the 
wavevector k is modified by rotating the sample around 
the vector q (Fig.fl^: the dipole matrix M is then given 
by M = {a) FR\a)^ where the matrix R represents a 

rotation of a degrees around the q axis. This expression 
corresponds to the predictions of Ref. m for an 5-wave 
modulation. 


Without loss of generality we choose to work in Carte¬ 
sian coordinates in which the incoming photon moves in 
the direction k = (0, 0,1) with polarizations a = (0,1, 0) 
and TT = (1,0,0). As shown in Fig. the polariza¬ 
tion of the outgoing photon are a' = (0, —1,0) and 
tt' = (cos 26>, 0, sin 20 ), and q = k' — k ~ (cos 6>, 0, sin 0 ). 
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For = 1, = 0 we then find : 


min/g/max/g ratios. 


( o ol cy, n \ ^ 

4 cos^ — sin — cos^ 0 sin Oj 

+ (cos^ a + sin^ a sin^ 0) ^ (A23) 

(a) = ^4 cos ^ sin^ ^ cos^ 0 sin 0 ^ 

+ (cos^ 0 — sin^ 6>(sin^ a + cos^ asin^ (A24) 


For Fz 0 one obtains more complex expressions, nu¬ 
merically depicted in Fig. If we normalize and 
by their maximal value (dotted curves of Fig. 14), 
we find that min/g/max/g ^ 0.75 for both 5 = 7r,cr. 
The experimental measurements (symbols) instead show 
min/g/max/g « 0.6. As explained by Comin et 
this discrepancy can be attributed to the self-absorption 
effect^. The corrected signal (dashed curves) gives a 
good agreement between the predicted and measured 


Comin et al additionally notice a significant asym¬ 
metry of the experimental data with respect to a ^ 
180° — a. According to a specific statistical model, they 
attribute this discrepancy to a dominant d-wave charac¬ 
ter of the charge modulation. Their classification implies 
the realization of a multi-band model with inequivalent 
Ox and Oy orbitals, in contrast to the present analysis 
based on a single-band model. Interestingly, the statisti¬ 
cal significance of their approach was not clearly estab¬ 
lished in the case of BSCCO. In addition. Fig, shows 
that the symmetry with respect to a ^ 180° — a is sig¬ 
nificantly modified by the corrections for self-absorption. 
We conjecture that the discrepancy between the present 
theoretical approach and the experimental measurements 
might be attributed to high-order corrections in the self¬ 
absorption, rather than to the symmetry of the charge 
modulations. 
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Note that current experiments would detect this situation 
as a pronounced peak at g = 2 /cf due to the common 
practice of “background subtraction”. 

The antinodes are here defined as the points at which the 
Fermi surface reaches the boundaries of the first Brillouin 
Zone 

Here we use the term “nesting” to indicate segments of the 
Fermi surface that are parallel to each-other, leading to 
an enhanced scattering at the corresponding wave-length 
difference 

The quantitative discrepancy between our predictions and 
their finding may be attributed to either the different ma¬ 
terial ( Lai. 88 Sro.i 2 Cu 04 ), or the different technique (res¬ 
onant elastic X-ray scattering). 

In the experimental curves reproduced in Fig. 06) > the 
signal in (0, q) and (g, q) directions were normalized by 


a different multiplicative factor Indeed, if the same nor¬ 
malization had been used, the two curves should had con¬ 
verged in the g ^ 0 limit. As a consequence, the plot of 
Fig. [^e) does not convey any information about the ra¬ 
tio between the (0.25,0) and (0.25,0.25) wavevectors, but 
only describes the behavior of the REXS intensity in the 
(0,g) and {q^q) direction, independently. 

Because the wavevector (^, 0) connects quasiparticles with 
opposite signs of the pairing gap, while the (g, q) wavevec¬ 
tor connects quasiparticles with the same sign, these two 
contributions are often referred to as “sign-reversing” and 
“sign-preserving”. 

Here we assume the impurity itself to be symmetric with 
respect to r —r. This assumption applies for example 
to impurities with s and d symmetry, but does not apply 
to impurities with p symmetry 


